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Pseudospin solitons in double-layer quantum Hall systems can be introduced by a magnetic field 
component coplanar with the electrons and can be pinned by applying voltages to external gates. We 
estimate the temperature below which depinning occurs predominantly via tunneling and calculate 
low-temperature depinning rates for realistic geometries. We discuss the local changes in charge 
and current densities and in spectral functions that can be used to detect solitons and observe their 
temporal evolution. 



73.40.Hm,71.45.-d,71.10.Pm,73.40.Gk,73.40.-c 

The study of multicomponent Quantum Hall sys- 
tems jjj has been enriched by the discovery of a variety 
of new phases. In double layers, the relevant discrete de- 
grees of freedom are labelled by the electron's layer and 
spin indices. At Landau level filling factor v — 2, recent 
theoretical work predicted several interesting phases ||] 
in which both layer and spin play a role; the existence 
of these phases has been confirmed by experiment jjj. 
The present study is on double-layer systems at filling- 
factor v — 1 0||. At this filling factor, the low-energy 
electron states are spin-polarized and the system has a 
broken symmetry ground state with spontaneous inter- 
layer phase coherence ||. The rich phase diagram for 
these systems, including the effects of in-plane fields, has 
been discussed at length in Ref. 

It is useful to describe this system using a pseudospin 
language j|,[3| in which pseudospin-up (-down) refers to 
an electron in the top (bottom) layer. The action is 
that of a two-dimensional ferromagnet with a hard-axis 
anisotropy and a Zeeman field perpendicular to the hard- 
axis The pseudospin configuration is specified by 
the spherical-coordinate fields 9(x,y), which describes 
the difference in charge density between the layers, and 
4>(x,y), which describes the relative phase of electrons 
in top and bottom layers. Phase solitons 4>o(x) exist as 
solutions to the classical equations of motion M. In this 
Letter we address the quantum dynamics of such soli- 
tons, predicting that, when pinned by applying gate volt- 
ages, depinning occurs at accessible temperatures pre- 
dominantly via quantum tunneling. This system offers 
a number of advantages for macroscopic quantum tun- 
neling studies, especially the possibility of using gate 
voltages and in-plane fields in combination to control 
metastable-state placement. We also discuss several lo- 
cal properties which can be used to detect solitons and 
observe their temporal evolution. 

Neglecting for the moment the effect of in-plane 
fields, the leading contributions to the imaginary-time 
effective action for the pseudospin field m(r,T) = 
(sin 9 cos (f>, sin 9 sin <j), cos 9), in the presence of both tun- 
neling and gates, is given by [|| 
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The first term is the Berry phase, conveniently expressed 
as rh • A = 0(1 — cos 6*). The gradients are the ex- 
change terms, with p the pseudospin stiffness. The term 
involving (3q is a hard-axis anisotropy, and in the follow- 
ing term, t is the tunneling amplitude which acts as an 
in-plane pseudofield. Finally, V(x) is the gate potential 
which can be adjusted in situ for appropriately fabricated 
double-layer samples. The local filling factor for each 
layer is V\ = (1 + cos6*)/2, v% = (1 — cos 9)/2, with the 
total filling factor v = V\ + i>i = 1 . The parameters p and 
t in Eq. ([!]) depend on m z . We require their expansion to 
quadratic order in i9 = 9 — tt/2, p(m z ) — po + p\-d 2 and 
t(m z ) = to + tid 2 , and shall use the Hartree-Fock results 
pi = -p and ii = -t /2. 

We limit our attention to finite-width (w) systems for 
which the soliton physics is particularly simple. The flex- 
ural bending mode JIH of a soliton in the x direction 
of a two-dimensional system has a spectrum given by 

e(ky) = 4:ir£ 2 k y ^J pq(2Pq + pok 2 ,), where k y is the trans- 
verse wavevector. The finite-size gap of these modes can 
be estimated by setting k y = n/w. For a given tempera- 
ture T, the quasi-one-dimensional limit results if flexural 
modes are frozen out, i.e., if the transverse sample size is 
smaller than § w(T) = {iz / t)[{ Pq / p ){l + ^/T+T 2 )] 1 l 2 , 
where r = fc^T '/ '(Att£ 2 (3q). For the sample parameters 
we give below, m;(1K) = 718 nm, whereas w(100mK) = 
7/Ltm. Assuming the transverse sample size to be less 
than this value, we can trivially integrate over the y di- 
rection. 

The estimates for the tunneling rates given below are 
based on the following (l(J parameter values: 2tt£ 2 Po w 
2.3 meV, p a w 0.024 meV, t w 0.1 meV, and^ = 11.8 nm. 
We limit our attention to the regime V(x) <C 2tt£ 2 (3q 
and look for solutions to the equations of motion with 
9 = tt/2. For V(x) = 0, the solution <j) (x - X) = 
4 arctan e^ x ~ x ^ s satisfies the equation of motion, where 
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S = y / 2irp /t £. The solution (j>o describes a static 2tt 
soliton with width 6 centered at position X. The pseu- 
dospins rotate in the xy plane by 2tt as we move through 
the soliton, beginning and ending at <f>o(±oo) = 0, 2tt. 
Quantum effects are incorporated by expanding about 
this classical solution, and dynamics is described through 
a canonical transformation to collective coordinates as 
we discuss below 0|. We include the effects of the 
gates self-consistently by substituting the soliton solu- 
tion (9,4>) = (7r/2,0o) into the full action. This is valid 
so long as the soliton energy £ 80 iiton = 8L y po/5 is much 
larger than the effective potential. 

Expanding to quadratic order in $(x, r) = 9(x,t) — 
7r/2, and performing the functional integral over $ (in 
the partition function), we obtain 
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where c = 2n£ y/8/3oPo is the spin- wave velocity. This 
equation is valid for t /(2ir£ 2 f3 ) <C 1 and Ap /(S 2 f3 ) -C 
1. For the parameter values given above, to/(2ir£ 2 /3 ) « 
0.04 and 4p /(6 2 f3 ) » 0.17. 

To describe motion of the soliton, we perform a canon- 
ical transformation to collective coordinates. This en- 
tails elevating X to a dynamical variable pj],p|, X — > 
X(t), and introducing a constraint to preserve the de- 
grees of freedom. We thus write 4>(x,t) — <po(x — 
X{j)) + ip(x — X(t),t) and expand to quadratic or- 
der in ip, where <po now describes a moving soliton and 
(p describes a dissipative spin-wave field, which is re- 
quired to be orthogonal to the zero mode of the soliton 
field. We incorporate this constraint, explicitly given by 
Jdx (j)' (x)ip(x,T) — 0, via the Fadeev- Popov technique. 
Briefly, the procedure || is to insert into the functional 
integral the identity fT>X S(Q) det(SQ/SX) = 1, where 
Q[X] — Jdxcj)' (x — X)<p(x,r). We then expand to sec- 
ond order in both tp and X/c. For the bounce solutions 
considered below, \X/c\ < 0.26. This yields a descrip- 
tion of the solitons, the spin waves, and the dynamic 
(non-linear) coupling between them. 

Integrating over the spin waves, and including the gate 
potential, yields an action corresponding to a particle 
(soliton) at position X, with bare mass M, moving in a 
potential V e g(X), with dissipation described by a nonlo- 
cal kernel K(t): 
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For low temperatures (fc^T <C Tiloq = 2hc/S), the kernel 
is given by K(t) = (1/tt) J °°duj J(w) e - W I T I. The spectral 



density J(u>) — (ui /loo8 2 )\Juj 2 — uJq Q(u> — luq) vanishes 
for u> < luq (8 is the step function), and the kernel K (r) 
is exponentially suppressed for |r| > l/^o- In the pa- 
rameter range we are interested in, fc^T <C huJb "C Tiwo 
(lui, is the bounce frequency — see below), we can expand 
X{t) - X(t') rs (t - t')X(t / ), and the effect of damp- 
ing reduces to a renormalization of the soliton mass M. 
This additive renormalization is intensive with respect 
to the transverse length L y of the sample whereas 
the bare mass is extensive (proportional to L y ), and so 
we shall not consider this additive piece in what fol- 
lows. For the estimates we give below, LOb/^o < 0.11, 
and kBTc/hu>b = 0.47, 0.16, 0.13, and 0.13 respectively 
for the entries listed in Table |, where Tc is the crossover 
temperature (see below). 

To leading order, the bare mass is given by M = 
h 2 Ly/(4Tr 2 £ 4 [3ofi)i and the effective potential by 
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This effective potential can be understood in terms of the 
change in the classical soliton energy density when the 
m z = values for t and p in Eq. (^) are replaced by their 
values at m z — —V(x)/4ir£ 2 f3o, i.e., when the pseudospin 
polarization produced by the bias potential is accounted 
for. In Eq. (§), we have neglected phase terms |J. All 
such terms vanish for incoherent tunneling, as in the case 
of tunneling out of a metastable state considered below. 

To obtain an explicit expression for V e s, we should 
specify the form of the gate potential V(x). We choose 
two simple square wells. One applies a gate voltage of 
—Vi over a width wi, —V2 over a width 102, and zero 
voltage elsewhere. The inside edges of the gates are a 
distance 2xq apart. This yields an effective potential ex- 
plicitly given by 



V cB (X) = v ^{-lyV? [tanhx, - tanh(^ 



i)}, (5) 



where v = -3L y p /[2(2Tr£ 2 ) 2 8p 2 }, z< = [X - (-l) l x ]/<5, 
and Wi = (—l) l Wi/6. Self-consistency demands that 
^soliton ^ V e fi, which places an upper bound on the gate 
voltage of ^ max w 2.3 meV. 

We will consider tunneling out of metastable states. 
Once a soliton tunnels out, it does not return. We 
therefore put \Vi\ < | V2 1 and W\ « i«2. The potential 
is schematically shown in Fig. [I], with the coordinates 
shifted as outlined below. The objective is to calculate 
the tunneling rate out of this metastable state. 

It is helpful to shift coordinates so that the metastable 
minimum is now at the origin and V c ff(0) = (See Fig. ^). 
There is a frequency ujq — y/V^(0)/M associated with 
the curvature of the metastable minimum, and a fre- 
quency LOb — y/V£fi(Xb)/M associated with the barrier. 
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FIG. 1. Effective potential produced by the gate voltages in 
the "metastable" configuration. The inset shows the "bare" 
gate potential. 

The tunneling rate is given by Y — Kloq (So/2-k) 1 ^ 2 e~ So , 
where if is a constant whose calculation requires explicit 
knowledge of the bounce trajectory. The exponent So is 
the action evaluated along the minimal (bounce) trajec- 
tory, which goes from X = to X = X e and back to 
X = 0, while r goes from — /3/2 to +f3/2, with f3 — > oo. 
The crossover temperature Tc defines the boundary be- 
tween thermally-dominated transitions (T > Tc) and 
tunncling-dominated transitions (T < Tc). Both transi- 
tions show exponential behavior and one usually assumes 
that the exponent dominates the prefactor, leading to the 
expression k B T c = V cS (X b )/So- 

The tunneling rate, along with the results for the 
crossover temperature Tc, the classical action So, and the 
attempt frequency ujq are listed in Table [l| for various gate 
voltage profiles and transverse sample sizes. These en- 
tries were evaluated using the typical model parameters 
quoted earlier which yield a soliton width of S = 14.5 nm 
and a bare soliton mass of M w m* (0AAL y /£), where m* 
is the conduction band effective electron mass in GaAs. 
We conclude from Table | that voltage profiles for which 
quantum tunneling will dominate at accessible tempera- 
tures are achievable with current submicron lithographic 
technology. 

There are at least three properties of pseudospin soli- 
tons which should make their existence and their motion 
experimentally observable: moving solitons disturb the 
charge balance between the layers; there exists a local 
electrical current circulating about the axis of the soli- 
ton; and the local quasiparticle gap is suppressed in the 
vicinity of the soliton. These local changes can be de- 
tected respectively by electrostatic, magnetic, and tun- 
neling probes of the 2DES. We discuss each briefly. 



TABLE I. Tunneling rate F, crossover temperature Tc, 
and attempt frequency ujq for several gate widths u>i, gate 
separations 2xo, and sample sizes L y . All entries have 
w 2 = 400 nm, Vi = -0.75 meV, and V 2 = -1.00 meV, ex- 
cept for the final entry, which also has V\ = — l.OOmeV. 
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That moving solitons transfer charge between layers 
may be seen by revisiting the classical equations of mo- 
tion for the action in Eq. (Q) and trying to find moving 
soliton solutions. Proceeding in powers of v/c, where 
v is the soliton velocity and c the spin wave velocity in 
Eq. (§), we write cl (s) = ir/2+{v/c)0 1 {s) + O(v 2 /c 2 ) and 
0ci(s) = <j>o{s) + (v I c)0i(s) + 0[v 2 1 (?) , where (j>o(s) is the 
previous solution for the static soliton, and s = x — vt 
is the coordinate of the moving frame. Note that this 
ansatz satisfies the equations of motion exactly for v = 0. 
To linear order in v/c, we find that <j>i(s) = 0, while 

±(c/2S)sech[( s -so)/6} 
01 [S> (a 2 A, + t ) - 6i sech 2 [(s - s )/6} ' l ° ' 

The local charge transfer produced by the passing soliton 
is related to the charge transfer produced by the static 
bias potential which, by virtue of the difference in layer 
energies it produces, gives the same rate of change of 
the interlayer phase 4>. This result shows that whereas a 
static soliton pseudospin rotates entirely in the xy plane, 
a moving soliton pseudospin rises out of the plane by an 
amount proportional to the velocity as it undergoes its 
2ir rotation. This is analogous to closely related models 
of anisotropic ferromagnets where exact moving soliton 
solutions are available p3[ . In physical terms, a moving 
soliton transfers charge between the layers, altering the 
electric potential profile outside of the sample, which can 
be measured using a scanning single-electron transistor 
[ ]l4| |. The total transferred charge for a soliton is given 
by Q e = L y ne J dx cos6* c i ~ ^irenLy (X/c)y/ 2po //?, with 
corrections of order {X/c) 2 and to/(a 2 f3o), and where 
n s» 10 n cm~ 2 is the electron density. Using the param- 
eter values given above, a system size of L y = 600 nm, 
and taking for X the maximal velocity along the bounce 
trajectory, we obtain \Q e \ ~ 2e. This value can be in- 
creased either by inducing a larger soliton velocity, or by 
going to a larger sample. These charges should appear, 
for example, as a soliton tunnels from a metastable state. 

Currents, much like those near a vortex in a super- 
conductor, circulate around the axis formed by a soli- 
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ton line and produce a magnetic field which could in 
principle be measured. The microscopic operator for 
the current flowing between layers is the time derivative 
of the layer-polarization operator in the Heisenberg rep- 
resentation. The only contribution to the commutator 
comes from the interlayer tunneling term in the micro- 
scopic Hamiltonian, and is proportional to the operator 
for the y component of the pseudospin. Taking its ex- 
pectation value we find for the 3D current between the 
layers (with the bottom to top layer direction taken as 
positive) jBT{x) — — t sin(0 o (x))/?i27r^ 2 . Currents flow 
in opposite directions on opposite sides of the soliton. 
The currents flowing within the bottom (B) and top 
(T) 2D layers are oppositely directed and given by ||] 
Jb{x) = —]t{x) = pd(po/h. For a stationary soliton 
these currents satisfy the local charge conservation con- 
dition djs(x)/dx = —djT{x)/dx = jBr(x). The current 
per unit length flowing around the soliton is ~ ep/hS 
and the characteristic loop area of the circulating cur- 
rents is ~ d5 where d is the layer separation. For typ- 
ical parameters, these currents produce a magnetic in- 
duction directed along the soliton axis with magnitude 
~ 0.1/xs/£ 2 (5, which could in principle be detected by a 
SQUID or with cantilever-based jl5| magnetometers. 

Microscopic calculations can also be used to evaluate 
the local gap for charged excitations as a function of posi- 
tion near a stationary soliton. We find pl| that the local 
gap is always reduced near the soliton center and is given 
by E gap (x) = £ , gap (oo) - 4t(l - cos(4> (x)). This change 
in spectral properties might be visible in scanning tun- 
neling probes of 2DES's. This gap suppression is local 
and will increase the local thermally activated quasipar- 
ticle resistivity p xx at the soliton center by a factor of 
exp(8t/ksT) ~ exp(10/T [Kelvin]) which can easily be 
very large. When a soliton tunnels away from a region 
defined by a set of voltage probes, we can expect a large 
signal to be seen in the conductances measured by those 
probes. 

In all of the above, we have assumed the external 
magnetic field was oriented perpendicular to the layers. 
But this need not be the case, and it may in fact be 
preferable to add an in-plane component to the field. 
In the presence of an in-plane field Bii = (0, — -Bii , 0), 
the energy [[[o| of the soliton is lowered and is given by 
E so \ = L y po(8/S - 47r 2 dB||/$o), where $ is the flux 
quantum. The benefit, in the present context, of an in- 
plane field is that it may be used to thermally create 
solitons at arbitrarily low-temperatures and thus ensure 
that solitons exist in the sample. 

In summary, we have shown that pseudospin solitons 
may be pinned and manipulated by applying external 
gates to the system. One can extend the analysis to the 
case of a double-well potential, which should exhibit an 
externally tunable tunnel splitting, and to the case of a 
periodic potential. The physics in the latter case should 
be driven by Berry and topological phases, and should be 



strongly affected — and possibly controlled — by the gate 
voltages. This opens the door to a variety of Aharonov- 
Bohm-type investigations, as well as soliton derealiza- 
tion and band formation ||. We have also shown that 
solitons have an electrical current circulating about their 
axis, that the local quasiparticle gap is suppressed in the 
vicinity of the soliton, and also that soliton motion in- 
duces a local charge imbalance between the layers. These 
properties may provide a means through which the soli- 
ton dynamics can be experimentally probed. 
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The expressions, correct to Hartree-Fock order, for the pseudospin stiffness p and tunneling amplitude t are given 
by p — Po and t = to- Consequently, the right-hand side of Eqs. (4) and (5) of our Letter should be multiplied by 1/2. 
Some of the entries in Table I are also affected and more representative values are given below. The conclusions and 
all other results remain unaffected. 

TABLE II. Tunneling rate F, crossover temperature Tc, and attempt frequency too for several gate widths wi, gate separations 
2xo, and sample sizes L y , assuming negligible dissipation (see text). All entries have W2 = 400 nm, Vi = — 0.75meV, and 
Vi — — l.OOmeV, except for the final entry, which also has Vi = — l.OOmeV (and the same Vi and W2 as the other entries). 
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